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ABSTRACT 

We present a new method to directly measure the opacity from H I Lyman limit (LL) absorption 
kll along quasar sightlines by the intergalactic medium (IGM). The approach analyzes the average 
( "stacked" ) spectrum of an ensemble of quasars at a common redshift to infer the mean free path A^jp 
to ionizing radiation. We apply this technique to 1800 quasars at z = 3.50 — 4.34 drawn from the Sloan 
Digital Sky Survey (SDSS), giving the most precise measurements on kll at any redshift. From z = 3.6 
to 4.3, the opacity increases steadily as expected and is well parameterized by A^p = Ao — b\(z — 3.6) 

with A = (48.4 ± 2.1)/if 2 1 Mpc and b x = (38.0 ± 5.3) Mpc (proper distance). The relatively 
high A^ f p values indicate that the incidence of systems which dominate kll evolves less strongly at 
z > 3 than that of the Lya forest. We infer a mean free path three times higher than some previous 
estimates, a result which has important implications for the photo-ionization rate derived from the 
emissivity of star forming galaxies and quasars. Finally, our analysis reveals a previously unreported, 
systematic bias in the SDSS quasar sample related to the survey's color targeting criteria. This bias 
potentially affects all z ~ 3 IGM studies using the SDSS database. 

Subject headings: large-scale structure of universe — quasars: absorption lines — intergalactic medium 



1. INTRODUCTION 

The observed high transmission of z ~ 3 quasars 
at rest wavelengths A r blueward of H I Lya reveals 
that the intergalactic m edium (IGM) is highly ionized 
(jGunn fc Peterson! 1 19651 ). The presence of the Lya for- 
est demands an intense, extragalactic ultraviolet back- 
ground (EUVB) radiation field. The quasars themselves 
provide a significant fraction of the required ionizing flux, 
buoyed by the emission from more numerous yet fainter 
star-forming galaxies. Several recent studies have ar- 
gued that t he latter population dominates the EUVB 
at z > 3 (iFaucher-Gieuere et alJ l2008al : ICowie et all 
l2009t iDairAglio et alJl2Q09D. where the quasar popula- 
tion likely declines (e.g. iFan et alJ 12004). These asser- 
tions, however, hinge on the opacity of the IGM to ioniz- 
ing radiation via H I Lyman limit absorption (kll) which 
directly impacts estimates of the EUVB measured from 
the integrated quasar and stellar ionizing emissivity. 

Traditionally, kll has been estimated from the in- 
cidence of so-called Lyman limit syst ems (LLSs) via 
surveys of quasar spectroscopy (e.g. lLanzettal 119911 ; 
IStorrie-Lombardi et alJ 11994 iPeroux et all 12001 ). Ob- 
servationally, one can rather easily identify systems with 
large optical depths T912 > 2 and the majority of these 
surveys have probed to this limit. The integrated opac- 
ity of the IGM, however, includes and is likely domi- 
nated by gas with rgi 2 < 1, the so-called partial Ly- 
man limit systems and Lya forest clouds. Systems with 
these H I column densities (iVgi ~ 10 14 — 10 17 cm~ 2 ) 
are especially difficult to survey because the strong 
lines of the Lyman series (e.g. Lya, Ly/3) lie on the 
flat portion of the curve-of-growth and they exhibit 
only weak absorption at the Lyman limit. Thcrc- 
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fore, current estimates of kll are based on an extrap- 
olation/interpolation of the frequency of systems with 
Nm< 10 14 cm- 2 and Nm > 10 1 7 5 cm- 2 (iMadau et ah 



-iVHT iu cm anu ivm j> iu cm i | iviauau et ai 
1999; ISchirber fc B ullock 2003; iFaucher-Giguere et al 



2009). Current constraints on «ll span over a magni- 



tude of uncertainty, especially at z > 4. 

In this Letter, we introduce a new technique to esti- 
mate kll that avoids the traditional line-counting statis- 
tics of the IGM. We analyze the average rest-frame spec- 
tra of 1800 z > 3.5 quasars drawn from the Sloan 
Digit al Sky Survey, Data Release 7 (|Abazaiian et al.l 
2009). As an ensemble at a common redshift, these 
"stacked" spectra show the exponential drop in flux at 
A < A912 = c/^gi2 = 911. 76A from the integrated opac- 
ity of the IGM. We precisely evaluate kll in a series 
of small redshift intervals covering z w 3.6 — 4.3 to ex- 
plore redshift evolution. We adopt a cosmology with 
H = 72/i 72 km s^Mpc -1 , Q m = 0.3, and ft A = 0.7 
and report proper lengths unless specified. 

2. METHODOLOGY 

The H I Lyman limit opacity of the IGM has tradition- 
ally been expressed as an effective optical depth r e s,LL 
estimated from the observationally constrained Nm fre- 
quency distribution f(Nm,z). An ionizing photon (y > 
t/912) emitted from a quasar with redshift z — z q will 
redshift to 1 Ryd at z = 2912 = (^912/^) (1 +Z q ) — 1- The 
effective optical depth that this pho ton experiences by 
IGM Lyman limit opacity is then (cf. iMeiksin &: Madaul 
1993): ' 



T e ff,LL(z912, Zq) 



(1) 



f(N m , z'){l - exp [ -N m a ph {z')]}dN m dz' 



2912 



where cr p h is the photoionization cross-section evaluated 
at the photon frequency. In practice, this approach is 
subject to large uncertainties because (1) f(Nm,z) is 
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Fig. 1. — The stacked spectrum (normalized at A r = 1450A) from 150 mock quasar spectra with IGM absorption derived from an 
assumed f(Nm,z) distribution. The adopted emission redshifts and S/N of the spectra were taken to match our SDSS quasar sample at 
z = [3.59,3.63]. The sharp decline in the flux at A r 912A is associated first with the higher order Lyman series opacity (for A r > 912A) 
and then Lyman limit continuum opacity (for A r < 912A). The dashed (cyan) curve overplotted on the spectrum is the predicted relative 
flux corresponding to T e g ll (Equation [TJ , normalized at A r = 912A. On top of this curve is the best-fit model of the flux evaluated over 
the interval A r = 848 — 905 A using our new methodology (Equation |6jl . The resulting estimates for the mean free path A^ f p are indicated 
on the figure and are in excellent agreement. Finally, the dotted curves show the predicted models after offsetting A^^ by ±10 Mpc and 
indicate the sensitivity of our technique to measuring A^J? . 



mfp ' 

poorly constrained for systems with T912 < 1, (2) ob- 
servational surveys rarely measure /(iVm, z) at the same 
epoch forcing interpolation and extrapolation; and (3) es- 
timates of /(-/Vhi, z) for the LLS may always be subjec t 
to large systematic uncertainty (jProchaska et al.l [2009) . 
This traditional approach may never yield a precise and 
robust estimate of r c ff LL- 

Our new approach is to directly measure T e s^ 
through analysis of averaged ensembles of quasar spec- 
tra. In Figure [TJ we present the stacked spectrum 
of 150 mock quasar spectra at z » 3.6. Each spec- 
trum was given a unique emission rcdshift z q and SED 3 
(normalized at 1450A), and then was blanketed with 
Lyman series and Lyman limit absorption from an 
assumed f(N m ,z) distribu tion (|Dall' Aglio et alj 120081 : 
IWorseck fc Prochas ka 2009) . The spectra were degraded 
to the nominal spectral resolution of the SDSS spectrom- 



eter (FWHM= 150km s 1 ) and Gaussian noise was added 
to give a distribution 4 of S/N values at A r = 1450A. The 
data was then averaged without weighting. 

Inspecting Figure [J one identifies the effective opacity 
of the Ly/3 forest at A r w 1000A and corresponding decre- 
ments in the spectrum at Ly7 and LyS. One then ob- 
serves a steep drop in the flux starting at A r » 920A due 
to the opacity of higher order Lyman series lines of op- 
tically thick absorbers (e.g. damped Lya systems). The 
continued decline at A r < A912, however, is dominated 
by the continuum opacity of H I. At all wavelengths, 
the scatter in the stacked spectrum is due to small-scale 
variance in IGM absorption, not noise in the individual 
spectra. 

Overplotted in Figure [1] is the flux model / = 
/gi2 exp[— T e s,hh(z)} with t^ll evaluated from the in- 
put f{Nui,z) distribution (Equation [JJ and /912, the 



3 Q uasar SEDs do no t show strong emission features at A r < 
912A dTelfer et al.|[2002T) . 



4 The S /N and emission redshift distributions correspond to the 
third bin in the observational analysis that follows. 
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flux at A = A912, estimated from the data. This is a 
good model of the stacked spectrum; even though the 
underlying average SED evolves as fx oc A 2 4 , the analy- 
sis is performed over too small a wavelength interval to 
note its evolution. The evaluation gives t^ll = 1 at 
(z q — 2912) = 0.22 corresponding to a proper mean free 



36.9 hj 2 L Mpc at z = 3.6. 



path A^ f 2 

Now consider an alternate evaluation of t^ll which 
follows the standard definition of optical depth, 



for a best-fit value n Zq = 0.028 hr^Mpc . The corre- 
sponding mean free path A^, = = 35.2 h? 2 Mpc 
is in excellent agreement with the traditional evaluation. 
We stress that the analysis was performed without any 
consideration of the quasar SEDs nor any consideration 
of evolution in the Lyman series line-opacity. Although 
these contribute to the observed flux in the stacked spec- 
trum, the exponential drop due to r e g ,ll dominates over 
this and any other astrophysical aspect. 



Teg,iz.(r,v) = / K LL (r',v)dr' 



(2) 



where the integral is evaluated to an arbitrary proper 
distance from the quasar. In an expanding universe, an 
ionizing photon emitted by the quasar will be attenuated 
by the Lyman limit opacity kll until it is redshifted to 
hv = 1 Ryd at z = Z912. If the photon is not absorbed 
by IGM line opacity from gas at z < Z912, it may be 
observed today at a wavelength A bs = (1 + 2912)0/1/912. 
During the photon's travel from z q to 2:912, the opacity 
«ll evolves because of the decreasing frequency (red- 
shift) and also from changes to the physical conditions 
of the universe (e.g. the expanding proper distance). We 
separate the frequency and radial dependencies in the 
opacity as follows, 



Khh(r, v) = K9ia(r) ( 

^912 



(3) 



where the frequency dependence related to <7 p h is ap- 
proximate. This treatment also ignores stimulated emis- 
sion, i.e. it assumes t^ll is dominated by 'clouds' with 
T912 < 1. 

In principle, one can adopt any radial dependence for 
kll- Expressing Equation [3] in rcdshift space, we have 



kll(z) = £912(2) 



l + z 
1 + 2912 



(4) 



With this functional form for the Lyman limit opacity, it 
is straightforward to integrate Equation [5] by adopting a 
Friedman- Walker cosmology where 



dr 

dz 



c/H 



H(z)(i + z) {i + z )^n m {i + z) 3 + n A ' 



(5) 



At z > 3, the universe is matter dominated and we can 

express dr/dz w c/(_ffo^m 2 )(l + z)~ 5 / 2 . Altogether, we 
find 

7"eff,LL(2912, Z q ) = 

z q 

^ T7? (l + 29i 2 ) 3 / k 912 (z')(l + zT 11/2 dz'. (6) 

Z912 

In practice, we find that the analysis of a single stacked 
spectrum does not constrain the redshift evolution in 
Therefore, we have simply parameterized kll by 
its value at z = z q , i.e., £912(2') = n Zq . The thin solid 
curve in Figure[T]shows the resulting flux model for T e ff ,ll 



3. RESULTS 

Our new approach provides a tight constraint on the 
effective mean free path near the quasar emission red- 
shift (at 2 w z q ). Because the stacked spectrum covers 
only several tenths in redshift below the Lyman limit, it 
imposes a very weak constraint on the redshift evolution 
of Instead, one must evaluate the stacked spectrum 
of quasars at a range of emission redshifts. 

We apply the methodology to 1800 quasar spectra 
drawn from the SDSS-DR7. We began with the vet- 
ted q uasar list from our s urvey of Lyman limit sys- 
tems ()Prochaska et al.ll2009f) which avoids all purported 
quasars in the SDSS-DR7 that have erroneous redshift 
estimates or are not bona-fide quasars. We have also 
ignored all quasars with strong associated absorption in 
the C IV, N V, and/or O VI doublets (e.g. broad ab- 
sorption line systems). We have not, however, removed 
quasars with evident Lyman limit absorption at 2 rj z q . 
It is our goal to estimate the entire Lyman limit opacity 
that quasars experience, except for the influence of gas 
on parsec scales. Therefore, our estimates of kll include 
opacity from the quasar's local galactic environment, i.e. 
its proximity region. 

The sample was further limited to the following cri- 
teria: (1) z q > 3.5 to insure significant coverage of the 
Lyman limit in the SDSS spectra and minimize the like- 
lihood that LLS bias the quasar target probability (but 
see below); (2) S/N > 4 at A r = 1450A; (3) z q < 4.35 
to insure that a stack of 150 quasars covers a redshift 
interval Az < 0.4. Starting at z = 3.5, we constructed 
a series of bins of 150 quasars each to produce a stacked 
spectrum. Each quasar spectrum was normalized by the 
observed flux at A r = 1450A (in a 20A window) and 
shifted to its rest-frame (nearest pixel). The full ensem- 
ble was then averaged (without weighting) ignoring bad 
pixels. A sample of three of the stacked spectra are given 
in Figure [2] 

As noted above, our analysis of r e ff ; LL includes con- 
tributions from the quasar's proximity region. In gen- 
eral, this corresponds to ~ 10 Mpc or AA r « 9A (e.g. 
iDall'Aglio et al.ll2008l ). From the stacked spectrum, Fig- 
ure [2] shows that the flux at the 'edge' of the proximity 
region (i.e. at A r « 900A) is w 0.75 times the flux at 
912A, giving Tefl^LL ~ 0.3. In the absence of the IGM 
beyond the proximity region, the flux would begin to re- 
cover shortward of 900A as v~ 3 . It is evident, therefore, 
that opacity from the 'true' IGM dominates our analy- 
sis. Furthermore, we find that our analysis yields models 
that extrapolate well to 912A suggesting the opacity of 
the proximity region follows the behavior of the general 
IGM. This is consistent with our finding that there is no 
strong differences in th e incidence of LLS near quasars 
(|Prochaska et al.ll2009h . 
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Fig. 2. — The stacked spectrum for three of our redshift bins (cyan: z = [3.59,3.63]; red: z = [3.86,3.92]; green: z = [4.13,4.34]) plotted 
against rest-frame wavelengths redshifted to the mean quasar redshift for the bin. Overplotted on each spectrum is our best-fit model for 
the absorbed flux below Agi2 due to Lyman limit opacity. These same curves are shown in the sub-panel against rest wavelengths. The 
emission lines are from Lyman series and metal-line transitions. 

Using the model for t^ll (Equation[5|), we have fitted 
the data to evaluate the opacity n Zq at z « z q in a series 
of redshift intervals, each containing 150 quasars. We 
minimized \ 2 over rest-wavelength intervals starting at 
905A (to minimize bias from strong Lyman series absorp- 
tion in the proximity zone of the quasar) and extending 
down in wavelength corresponding to the larger of 3900A 
and z q — zgi2 = 0.4. The scatter in the stacked spectrum 
was estimated locally in 21 pixel bins centered at each 
data point and presumed to be Gaussian. The flux at 
A912 was estimated from the data but allowed to vary by 
10% when minimizing y 2 . Because the observed scatter 
in the stacked spectrum is systematic (related to stochas- 
ticity in the IGM) and is not included in our model, one 
cannot estimate a(n Zq ) the uncertainty in n Zq from stan- 
dard x 2 techniques. Instead, we performed a bootstrap 
analysis of 100 realizations of each stacked spectrum and 
estimated <j(n Zq ) from the resulting distribution of k z 
values. 

Figure [3] presents the evaluations of K Zq in terms of the 
mean free path (AjjJ^, = 1/ K Zq ) . The mean free path ex- 
hibits a peak (minimum in opacity) of nearly 50 fi^ 2 Mpc 
at z — 3.6 and declines with increasing redshift. We 



have parameterized the redshift evolution in A^, with 
a simple linear regression: X^ p — Xq — b\{z q — 3.6). 
Restricting the analysis to z > 3.59, a y 2 minimization 
of this model to the binned evaluations of A^, gives 

A = (48.4±2.1)/i 72 1 Mpc and b x = (38.0±5.3) /i^Mpc. 

Formally, our analysis also indicated a rise in the opac- 
ity at z < 3.6, i.e. the lower limits in Figure [3J This runs 
contrary to all expectation and current understanding of 
the IGM. Initially, we suspected that this measurement 
indicated a systematic e rror in the SPSS spec tra at the 
bluest wavelengths (e.g. iBernardi et al.l [2003) . To test 
this hypothesis, we examined the (u — g) colors of the 
quasars in the first (z w 3.5) and third (z sa 3.6) quasar 
bins. Figure [?] histograms the two (u — g) distributions. 
The colors of the z = 3.5 quasars are systematically red- 
der than those at z = 3.6; this is the exact opposite 
of what one predicts if the IGM were monotonically in- 
creasing in opacity with redshift. We conclude that the 
quasars at z = 3.5 drawn from the SDSS are redder than 
the cohort at z = 3.6 because of an elevated incidence of 
Lyman limit opacity, confirmed by a nalysis on the inci- 
dence of LLS (|Prochaska et al.ll2009D . 

We then explored whether this elevated opacity is re- 
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FlG. 3. — Estimates of the mean free path (l/K Zq ) versus redshift. Error bars were estimated from standard bootstrap analysis. 

The red curve is a simple linear regression of the binned data: = A - b x (z - 3.6) with A = (48.4 ± 2.1)hf 2 1 Mpc and b x = 

(38.0 ± 5.3) h~2 Mpc. These are the most precise estimates of A^j? to date and provide the first robust description of its redshift evolution. 
The values at z < 3.6 (shaded region) are severely affected by a systematic bias in the SDSS database (see text) and are presented here 
only as lower limits. 



lated to observational bias in the SDSS quasar sam- 
ple. We simulated the SDSS experiment by construct- 
ing mock quasar spectra at z ks 3.5 and z « 3.6 with 
intrinsic SEDs having mean (u — g) color of 0.57 mag 
and standard deviation of 0.19 mag. These spectra were 
blanketed with IGM absorption assuming a monotoni- 
cally increasing opacity with redshift. After restricting 
the quas ar sample according to the SDSS color-selection 
criteria ([Richards et al.l 12002ft . we found the z ss 3.5 co- 
hort has systematically redder colors than the higher red- 
shift sample. Similarly, we find a correspondingly higher 
opacity inferred from the stacked spectrum. Because of 
the targeting criteria, the cohort of z « 3 — 3.6 quasars 
in the SDSS spectroscopic database are systematically 
biased against having (u — g) < 1.5 which biases against 
sightlines without strong Lyman limit absorption. Our 
analysis indicates that the bias extends to z q = 3.6, be- 
yond which very fe w quasars are predicted to e ver have 
such blue color (see iWorseck fc Proch aska 200j| for fur- 
ther details). The results for z > 3.6 are presented in 
Table □ 

4. DISCUSSION 



TABLE 1 
SUMMARY TABLE 



z 


<z q > 


•^analysis 


(Mpc^ 1 ) 


CT ( K z g ) 






(A) 


(Mpc" 1 ) 


[3.59,3.63] 


3.61 


[846, 905] 
[839, 905] 


0.0209 


0.0023 


[3.63,3.67] 


3.65 


0.0218 


0.0021 


[3.67,3.71] 


3.68 


[834, 905] 


0.0224 


0.0026 


[3.71,3.76] 


3.73 


[835, 905] 


0.0212 


0.0020 


[3.76,3.81] 


3.78 


[836, 905] 


0.0257 


0.0020 


[3.81,3.86] 


3.83 


[837, 905] 


0.0254 


0.0024 


[3.86,3.92] 


3.89 


[837, 905] 


0.0243 


0.0024 


[3.92,4.02] 


3.96 


[839, 905] 


0.0292 


0.0026 


[4.02,4.13] 


4.08 


[840, 905] 
[842, 905] 


0.0341 


0.0033 


[4.13,4.34] 


4.23 


0.0403 


0.0036 



Note. — Because of the color-criteria bias discussed in the 
paper, we caution that the values for the first few bins may sys- 
tematically underestimate K Zq by 10 — 30%. All wavelengths are in 
the quasar rest-frame and all distances are proper. The assumed 
cosmology has Q m = 0.7, Q\ = 0.7, and H$ = 72 km s- 1 Mpc — 1 . 

Our new technique provides the most precise mea- 
surements on the IGM opacity to H I ioniz- 
ing radiation at any redshift. Previous estimates 




Fig. 4. — Histograms of the (u — g) colors of quasars in the redshift intervals z = [3.50,3.55] (black) and z = [3.59,3.63] (red). The 
lower redshift cohort shows systematically redder (u — g) colors and a two-sided Kolmogorov-Smirnov test rules out the null hypothesis 
at greater than 99% c.l. This offset in (u — g) color contradicts standard expectation of a universe where H I absorption monotonically 
incre ases with redshift. It occurs because of a bias in the SDSS quasar sample related to the survey's targeting criteria lIRichards et al.l 
l200l) . Note, objects with (u - g) > 4.5 are generally not detected in the u-band. 



of kll were limited by large uncertainties in the 
iVni frequency distribution, especially t he incidence 

of systems with Nm ~ 10 17 cm~ 2 (iMadau et al.l 

19991 ISchirber fc Bullockl 12001 [Mhalda-Escudel Hool 
Meiksin fc White) 12004 IFaucher-Giguere et alj l2008al ). 
Our results indicate that systems dominating kll at z > 
3.6 must evolve more slowly than the TVhi < 10 14 cm~ 2 
Lya forest. In turn, we infer a flattening in the H I 
frequency distribution between TVhi ~ 10 15 — 10 17 cm -2 
(jPetitiean et alj 119931 : lO'Meara et~ai1 I2007D . We pro- 
vi de a full discussion on the implications for f(Nm,z) 
in lProchaska et all (j2009l ). 

Our analysis also gives the first direct description of 
the evolution in A;^ f 2 , albeit over a small redshift inter- 
val. Our results are well parameterized by a linear de- 
crease in A^ f p but can also be described by a (1 + z)~ 7 
power-law with 7 = 3.5 — 5.5. These values are consis- 
tent with the observed evolution in the incidence of LLS 
(jProchaska et al.ll2009l ). 

Our results refine re c ent inferences (e.g . 
IFaucher-Giguere et al.l l2008al : lDaH'Aglio etTal! 120091) 
that galaxies contribute significantly to the EUVB at 
z > 3. We assume the photoionization rate at z = 4 



inferred from the effec tive Lya opacity of the IGM 
(logTiGM = -12.3; IFaucher-Giguere et alj l2008bh . 
Comparing this value against the photoionization 
rate inferred from quasars T q using an emissiv- 
ity Cq 19 = 2 x 1 24 ergs~ 1 Hz~ 1 Mpc~ 3 (comoving; 
iHopkins et all l2007t iBongiorno et alj 120071) and adopt- 
ing our estimate of A^p at z — 4, we find T q = 0.5 Tigm- 
This suggests a modest but non-negligible contribution 
from galaxies to the EUVB at this redshift. The sys- 
tematic uncertainties in £g 12 and Tiqm are sufficiently 
large that one could recover riQM=rq, but T q would 
overpredict Tigm at z — 3.5 given the observed rise in 
^mfp an d e 9i2 w ith decreasing redshift. 

The results also revise at least s ome p revious estimates 
of the EUVB. lHaardt fc Madaul (|1996f) and their subse- 
quent analyses (CUBA), for example, have adopted an 
approximately 3 times shorter mean free path at z ~ 4 
than our analysis reveals. This implies: (1) the nor- 
malization of their EUVB spectrum is several times too 
low; and (2) the EUVB spectrum is softer at energies of 
w l,Ryd. The latter point may help to reconcile appar- 
ent contradictions in the metal-line analysis of the IGM 
(jAguirre et al.l 120011 ) without resorting to a large input 
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from galaxies. Other analyses, however, have used esti- 
mates for A^fp that are in much better agreement with 
our results (e.g. lFaucher-Giguere et al1l2009f h An under- 
standing of the full implications of our new constraints 
on A^ f p awaits new calculations of the EUVB. 

We identified a previously unreported systematic 
bias in the SDSS quasar spectroscopic sample against 
sightlines at z < 3.6 that are 'clear' of optically 
thick absorbers. This bias has implications for a 
range of IGM analysis at z ~ 3 including: (i) the 
paucity of sightlines for st udying He II reionization 
(|Worseck fe Pr ochaska 2009 ); (ii) an overe stimate of the 
incidence of LLS dProchaska et alj I2009D and damped 
Lya systems ([Prochaska fc Wolfe! I2009D ; and (iii) analy- 
sis of the Lyct forest. We caution that all existing studies 
of the IGM at z ~ 3 using the SDSS database should be 
reviewed in light of this systematic bias. 

In future work, we analyze these same spectra to con- 



strain the SEDs of z > 3.5 quasars, infer the Nm fre- 
quency distribution and Doppler parameters of gas with 
Nai w 10 16 cm~ 2 , and isolate the opacity of the IGM 
far from the quasar's proximity region. The technique 
introduced here is easily extended to higher and lower 
redshifts by obtaining modest signal-to-noise (S/N), low- 
resolution spectroscopy of several hundred quasars. Fu- 
ture ground and space-based programs will precisely es- 
timate kll from z q rs 0.5 — 5. 
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